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ABSTRACT 

The recently approved Juno mission will orbit Jupiter for one year in a highly 
eccentric (rmin = 1-06-Rjup, Tmax = 39i?jup) polar orbit {i = 90 deg) to accu- 
rately map, among other things, the jovian magnetic and gravitational fields. 
Such an orbital configuration yields an ideal situation, in principle, to attempt a 
measurement of the general relativistic Lense-Thirring effect through the Juno's 
node Q which would be displaced by about 570 m over the mission's duration. 
Conversely, by assuming the validity of general relativity, the proposed test can 
be viewed as a direct, dynamical measurement of the Jupiter's angular momen- 
tum S which would give important information concerning the internal structure 
and formation of the giant planet. The long-period orbital perturbations due to 
the zonal harmonic coefficients Je, £ = 2, 3, 4, 6 of the multipolar expansion of the 
jovian gravitational potential accounting for its departures from spherical symme- 
try are, in principle, a major source of systematic bias. While the Lense-Thirring 
node rate is independent of the inclination i, the node zonal perturbations vanish 
for i = 90. In reality, the orbit injection errors will induce departures 6i from the 
ideal polar geometry, so that, according to a conservative analytical analysis, the 
zonal perturbations may come into play at an unacceptably high level, in spite 
of the expected improvements in the low-degree zonals by Juno. A linear combi- 
nation of Q, the periJove u and the mean anomaly Ai cancels out the impact of 
J2 and Jq. a two orders of magnitude improvement in the uncanceled J3 and J4 
would be needed to reduce their bias on the relativistic signal to the percent level; 
it does not seem unrealistic because the expected level of improvement in such 
zonals is three orders of magnitude. More favorable conclusions are obtained by 
looking at single Doppler range-rate measurements taken around the closest ap- 
proaches to Jupiter; numerical simulations of the classical and gravito-magnetic 
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signals for this kind of observable show that a 0.2 — 5% accuracy would be a 
realistic goal. 

Subject headings: Relativity and gravitation; Lunar, planetary, and deep-space probes; 
Jupiter; Gravitational fields 
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Introduction 



20071 ) to Jupiter. Juno 



Recently, NASA approved the Juncu mission ( iMatousekl 
is a spinning, solar powered spacecraft to be placed in a highly eccentric polar orbit 
around Jupiter (see Table [T] for its relevant orbital parameters) specifically designed to 
avoid its highest radiation regions. Understanding the formation, evolution and structure 
of Jupiter is the primary science goal of Juno. It will carry onboard a dual frequency 

Table 1: Planetocentric nominal orbital parameters of Juno, a, e, i are the semi- major axis (in jovian radii 
R = 71492 km), the eccentricity and the inclination (in deg) to the Jupiter's equator, respectively. P is the 
orbital period (in days). T is the mission duration (in years). 



a (R) 



z(deg) P(d) T(yr) 



20.03 0.947 90 



11 



gravity/radio science system, a six wavelength microwave radiometer for atmospheric 
sounding and composition, a dual-technique magnetometer, plasma detectors, energetic 
particle detectors, a radio/plasma wave experiment, and an ultraviolet imager/spectrometer. 
The nominal mission's lifetime is 1 year. J uno is aimed, am ong other things, at accurately 



mapping the grayitati onal field of Jupiter (lAnderson 



([Anderson et al. 



19761 ) with unprecedented accuracy 



20041 ) by exploiting the slow apsidal precession of its 11-day orbit. 



In this paper we wish to explore the possibility offered by Juno to perform a test of 
general relativity by directly measuring the gravito- magnetic Lense-Thirring effect; its 
basics are reviewed below. Even putting aside the more or less successful attempts so far 
performed with other natural and artificial test particles orbiting different central bodies, it 
must be recalled that a satisfactorily empirical corroboration of a fundamental theory like 



^See on the WEB 



http://juno.wisc.edu/index.html 
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general relativity requires that as many independent experiments as possible are conducted 
by different scientists in different laboratories; thus, it is worthwhile to try to use different 
gravitational fields to perform such a test of intrinsic gravito-magnetism. Conversely, since, 
as we will see, the Lense-Thirring precessions are due to the proper angular momentum S 
of the orbited central body, one may also assume the existence of the general relativistic 
gravito-magnetism and consider such a test as a direct, dynamical measurement of the 
Jupiter's angular momentum through the Lense-Thirring effect; this would yield further, 
important information concerning the interior of Jupiter. Indeed, the moment of inertia 
ratio CiM lf' entering 5* is a measure of the concentration of mass towards the center of 



the planet (jlrwin 



2003 ). Such a figure, together with the measured values of the zona^ 



coefficients of the gravity field accounting for its deviations from spherical symmetry may 



composition vary with depth (llrwin 



Irwin 


2003; 


Guillot 


2005; 


Hori et al. 


2008) 



dynamical, model-independent determination of S would be important also for a better 



knowledge of the history and formation of Jupiter ( iMachida et al. 



20081). 



L et us, now, briefly review the basics of the Lense-Thirring effect (ILense and Thirring 



1918 



Zel'dovich and Novikov 



1971 



Soffel 



19891 ). In its weak-field and slow- motion 



approximation, the field equations of general relativity get linearized looking like those 
of the Maxwellian electromagnetism. Analogously with the magnetic field generated by 



moving elec 



(IMashhoon 



;ric c harges, mass-energy currents give rise to a gravito-magnetic field 

ar from an isolated spinning body of mass M and proper angular 



20071 ) 



momentum S it is (ILichtenegger and lorio 



2003) 



B. = --^ [5-3(5. f)f]. 



cr- 



where G is the Newtonian gravitational constant and c is the speed of light in vacuum. B^ 



^They preserve the axial symmetry. 
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exerts the non-central Lorentz-like acceleration (jSoffellll989l : iMashhoonl 120071 ) 



iGM 



— V X B„ 



(2) 



upon a test particle moving with velocity v. For ordinary astronomical bodies like, 
e.g., a planet, A'^^ is many orders of magnitude smaller than the Newtonian monopole 

= GM/r"^, so that it can be considered as a small perturbation. As a result, the 
longitude of the ascending node f l and the argument of pe ricentre a; of a test particle 



undergo small secular precessions (ILense and Thirring 

2GS 



19181) 



LT 



(3) 



6GS cosi 



CULT 



g2^3/2 ■ 



(4) 



Concerning a direct measurement of the Lense-Thirring effect, attempts have been 



recently performed wit 



(ICugusi and Proverb 



I the LAGEOS satellites in the gravitational fie ld of the Earth 



1978 



Ciufolini and Pavlis 



Surveyor probe or biting Mar s ( 



the Solar System 



lorio 



2006 



2004 



Riesetal 



2008f ). the Mars Global 



lorio 



(jCiufolini and Pavlis 



2005 



2009al ) and some of the inner planets of 



2008I). The e valuation of the total accuracy of th e LAGEOS 



lorio 



2007 



lorio 



2009al ) tests has raised 



2007al ) and MGS flKroghl l: 

a debate in the recent past because of the difficulty of realistically assessing the impact of 
certain competing dynamical effects acting as sources of systematic errors; f or example, the 



total accuracy of the LAGEOS test may be as large as some tens p ercent ( 



similar shortcomings may affect also the approved LARES mission (jlorio 



lorio 



2009c 



2009bh : 



d|). For 



an ov erview of se veral theoretical and experimental features of gravito-magne tism see 



e.g., (jlorio 



2007bl ). Concerning the jovian scenario. 



Lense and Thirring ( 119181 ) originally 



proposed to use the orbital p recessions of the Galile an satellites; such a possibility has 



been recently i n yestig ated by 



lorio and Laineyl ( 120051 ). but it seems to be still premature. 



Haas and Rossi (119751 ) proposed a spacecraft-based experiment in the gravitational field of 
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Jupiter to measur e another gravito-magnetic effect, i.e. the precessio n of a gyrosco pe (jPugh 



1959 



Schiff 



1974 ) whose 



19601 ). It is also one of the goals of the GP-B mission f lEveritt 
target was a ~ 1% measurement of such an effect with four superconducting gyroscopes 
carried onboard by a low-altitude polar spacecraft in the gravitational field of the Earth, 
but it is still unclear if it will be finally possible to meet the original accur acy because of 



some unexpected systern atic aliasing effects occurred during the missioio ( IConklin et al. 



2008 



Everitt et al. 



20091 ). A test of gravito-magnetisno concerning the deflection of 



electromagnetic waves by Jupiter in its orbital motion has been performed in a dedicated 



radio-interferometric experiment ( iFomalont and Kopeikin 



20081 ). With regard to other 



suggested non-gravito-magne tic tests of general relativity in the jovian gravitational field. 



Hiscock and Lindblom 



pericenter precessions ( lEinstein 



19791) propose d to measure the much larger gravito-electric Einstein 



19151 ) of the natural satellites of Jupiter and Saturn. 
There exist also plans for performing a test of the light bending due to the Jupiter's 
monopole and quadrupole mass moments with the forthcoming astrometric mission GAIA 



( ICrosta and Mignard 



20061). 



The Jupiter's proper angular momentum amounts to (jSoffel et al. 



20031 ) 



5 ^ 6.9 X 10^^ kg m^ s~\ 



(5) 



^See on the WEB http://einstein.stanford.edu/ 

'^In this case, the mass currents inducing a gravito-magnetic action are not those related 
to the Jupiter's proper rotation (intrinsic gravito-magnetism) , but are due to its translational 
orbital motion (extrinsic gravito-magnetism). 
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Table [T] and eq. ([5]) yield for Juno 

^LT = 68.5 mas yr~^ (6) 

COLT = 0. (7) 
which correspond to a shift Au of the cross-track component of the planetocentric position 



(jChristodoulidis et al. 



19881) 



Az/LT = a-y/ 1 + ^ sin zA^lt = 572 m (T = 1 yr) (8) 

over the entire duration of the mission. A total accuracy of the order of 1-10 m with 
respect to the km-level of the past Jupiter missions in reconstructing the Juno's orbit in 
a planetocentric frame does not seem an unrealistic target, although much work is clearly 
required in order to have a more firm answer. Note that a 1-10 m accuracy implies a 
0.2 — 2% error in measuring the gravito-magnetic shift 

The fact that the possibility of detecting the Lense-Thirring effect with Juno's orbit 
seems worth of further consideration can be preliminarily shown also with a different 
approach with respect to the cumulative measurement over the full mission duration 
previously outlined. Indeed, a gravity-science pass for Juno is defined by a continuous, 
coherent Doppler range-rate measurement plus and minus three hours of closest approach; 
in practice, most of the Lense-Thirring precession takes place just during such a six-hours 
pass, a near optimum condition. Another crucial factor is the orientation of the Earth 
to the Juno's orbit: our planet will be aligned 67 deg from the probe's orbital plane at 
approximately two degrees south latitude on the jovian equator. Preliminary numerical 
simulations of the Juno's Lense-Thirring Doppler range-rate signal show that such an 
orbital geometry represent a perfect compromise for measuring both the Jupiter's even zonal 
harmonics and the gravitomagnetic signal itself. Indeed, it turns out that the maximum 
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Lense-Thirring Doppler signal over a single six-hour gravity pass is of the order of hundred 
fim s~^, while the limit of accuracy for Juno's Ka-band Doppler system is about one fim 
over such a pass. Thus, even by taking 25 repeated passes, over a total of approximately 
33, it would be possible to reach a measurement precision below the percent level. 

Concerning the magnitude of the jovian gravito-magnetic field, it must be noted that 
in literature there are estima tes for 5* which point towards smaller values than eq. (|5]) by a 



factor 1.5 — 1.6; for example. 



Machida et al. 



fl2008[ ) yield 



S = 4.14 X 10^^ kg m^ s"^; 



(9) 



the ratio of eq. (|5]) to eq. 



i.e. 1.6, is close to 1.5 coming from the ratio o: 



C/MR^ = 2/5 = 0.4, vahd for a homogenous sphere, to C/MR^ = 0.264 bv llrwinl ( 120031 ) 
who assumes a concentration of mass towards the Jupiter's center. Here we will consider 
only the systematic uncertainty induced by the imperfect knowledge of the Newtonian part 
of the Jupiter's gravitational field; we will use eq. for S. 



2. The systematic uncertainty due to the zonal harmonics of the Jupiter's 

gravitational potential 



A major source of systematic uncertainty is represented, in princip 



of the Jupiter's gravitational field from spherical symmetry (lAnderson 



e, by the departures 



19761). 



2.1. Analytical calculations 

Indeed, the zonal (m = 0) harmonic coefficients Ji of the multipolar expansion of 
the Newtonian part of the planet's gravitational potential give rise to long-period, that is 
averaged over one orbital revolution, orbital effects on the longitude of the ascending node 
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Q, the argument of pericentre u and the mean anomaly Ai ( lKaulalll966l ) 



(10) 



i=2 



1=2 



M 



i=2 



(12) 
(13) 



where VL,i,Cj,i,M.,i are coefficients depending on the planet's GM and equatorial radius 



i?, and on the spacecraft's inclination i and ecc entricity e t 



and eccentricity Gipq{e) functions, respectively (jKaulalll966l ). Note that one of the major 



irough the inclination Fi^npii) 



scientific goals of the Juno mission is a greatly improved determination of just the harmonic 
coefficients of the jovian gravity potent ial; for the prese i it-day values of the zonalqj for 



2,3,4,6 see Table [21 According to 



Anderson et al. 



(120041 ). it might be possible to 



Tab 



2003!) 



e 2: Zonal harmonics of the Jupiter's gravity field according to the JUP230 orbit solution (jjacobson 



Js (xlO^ 



J3 (xlO^) J4 (xlO^) Je (xlO^) 



14696.43 ± 0.21 -0.64 ± 0.90 -587.14 ± 1.68 34.25 ± 5.22 



determine the ffist three even zonals with an accuracy of 10 ^ and the other ones up to 
£ = 30 at a 10"^ level. Concerning J2, this would be an improvement of two orders of 



^The Jupiter gravity field is essentially determined by the Pioneer 11 flyby at 1.6-Rjup 
( lAnderson Ill976l ): Voyager added little, and Galileo, which never got close to Jupiter, added 
nothing. 
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magnitude with respect to Table [2l while the improvements in J4 and Je would be of the 
order of three orders of magnitude. By using the results we are going to present below for 
the long-period node and pericenter precessions, it can be shown that determining the low 
degree zonals at 10~^ level of accuracy translates into an accuracy of the order of 0.5 — 1 
mas yr~^ in Q and u, thus confirming the expectations of the previous Section. 

For the long-period terms the condition 



2p + q = 



(14) 



is fulfilled. Thus, 



n 



R 



\J\ — sin i\a- 



E 

p=0 



dFta 
di 



'P 



(15) 



n 



R 



p=0 



cot i- 



dF, 



di 



G 



ep{2p-£) + 



Ffi 



dC 



ip(2p-i) 

de 



W^^0p{2p-^),(16) 



M 



n < 1 



p=0 



6G. 



lp{2p-£) 



(1 



dG 



ep(2p-e) 



de 



eop{2p~i) 



(17) 



where WeQp(^2p-e) are trigonometric functions having the pericentre as their argument and 
n = ^GM/a^ is the unperturbed Keplerian mean motion. Contrary to the case of small 
eccentricity satellites, in this case we will be forced to keep all the terms of order 0{e^) 
with A; > 2 in computing the eccentricity functions for given pairs of £ and p. Moreover, we 
will need also all the non-zero eccentricity and inclination functions for a given degree i, 
that is we will consider all the non- vanishing terms with < p < i. First, we will extend 



our calculations to the even zonals so far determined, i.e. 



2, 4, 6. In this case. 



Wiop(2p-i) = cos[{£ - 2p)ujt] = cos(gcjt). 



;i8) 



It must be noted that the terms with 



P 



2.« = 



(19) 



yield secular precessions, i.e. W^q-o ~ while those with 

q = 2p-£^0 (20) 

induces harmonic signals with circular frequencies —qu. 

For the degree i = 2 the non- vanishing inclination and eccentricity functions and their 
derivatives are 

F2oi = ^sin2^-i. (21) 



dF2oi 3 . . . 

— ; — = -sinzcosz. (22) 
di 2 ^ ^ 



*^2io - _g2)3/2- (23) 



dG2io 2 



de 3(l-e2)4/3- 
In this case, only secular precessions occur. 
For £ = 4 we have 



(24) 



35 15 

^401 = -— sin^z + — sin^z = F403, (25) 
62 lo 



105 4 15 2 3 
^402 = — sm 2 - y sm 2 + -. (26) 



dF^oi 35 . 3 15 . A . dF^os 

— ; — = sm z H smz cosz = — ; — , [27) 

di V 8 8 / di 



dF,02 /105 . 3. 15 . . 

— ; — = sm I -smz COS2. [2t 

di V 16 4 ' 
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3 

^^41-2 = 4^^_g2)7/2 = ^432, (29) 



1 + fe^ 

dGii-2 _ 3 e (7 - 6e^) _ 6^^432 \ 
~ 14(l-e2)8/7 ~ rfe ' ^ 

dG42o 1 e (23 - ISe^) 



de 7 (l-e2)8/7 

In addition to secular terms, also harmonic signals with the frequencies ±2a; are present. 

In the case of £ = 6 the inclination and eccentricity functions, along with their 
derivatives, are 

693 . 315 



(32) 



601 



602 



603 



sin^ i — sin"^ i = -Fgns, (33) 
512 256 ' ^ ' 



3465 . 6 . 315 . 4 . 105 . 2 . ^ 

sm I H — sm i — sm i = reoA^ (34) 



1024 64 64 



1155 . 6 . 945 . 4 . 105 . 2 . 5 ,^ , 

sm t sm t -\ sm t . 35) 

256 64 32 16 ^ ^ 



dFeoi ^2079 ^ 5 . 315 . 3 A . _ dF^o^ 



sm t — 7-r sm ? I cos? = — (36) 



di V 256 64 y di ' 



dFeo2 f 10395 . 5 . 315 . 3 . 105 . \ . ^^504 

sm « H sm z smz cosz = — — — , (37) 



di V 512 16 32 / di 



dFeo3 1^3465 . 5 . 945 . 3 . , 105 



di V 128 16 16 



sin i sin i H sin i cosi. (3^ 
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5 

^61-4 = Y^7Y3^5jTT7^ = <^654, (39) 



5 / e^^ 

^62-2 = TTTq QVnT? (■'-"'" "o" ~ ^"642, (40) 



2 (1 - e2)iV2 \^ 8 



G630 = T. + + ^eM . (41) 



;i - e2)ii/2 



dGQi-4 5 (l — ||e^) dG^^^ 



de 4 (1 - e2)i2/ii c^e 



(42) 



de (i_e2)i2/ii ' ^ ^ 

dGeso ^ 7 e (64 - lOe^ - 45e^) 

44 (l-e2)i2/ii ' ^ ^ 

In addition to the secular rates, also harmonic signals with frequencies ±4a;, ±2w are 
present. 

Let us, now, focus on the action of the odd {i = 3, 5, 7, ...) zonal (m = 0) harmonics. 
In this case, 

Wa)p^2p-e) = sm[{i - 2p)ujt] = - sin(gwt), (45) 
so that only harmonic terms occur for q ^ 0. 



For £ = 3 we have 



-IS- 



IS 3 

^301 = —sm^i--smi = -F302, (46) 
ib 4 



dFsoi /4S . 2 . 3\ . dF302 

— ; — = — sm t cosz = ; — . (47) 

di \ 16 4 di ^ ^ 



- _ g2)5/2 - ^^321, (48) 



de (l-e2)6/5 de ' 

Thus, we have long-period effects varying with the frequencies ±u. 



(49) 



Concerning the even and odd zonal long-period harmonic terms, it must be noted that 
they can be approximated by secular precessions with a high level of accuracy over the 
expected 1-yr lifetime T of Juno because the period of its periJove is of the order of ~ SOO 
yr, i.e. 

cos{quj) = cos{quJt + quq) ~ cos(ga;o), < t < T, (SO) 

sin a; = sin(a;t -|- uq) ~ sincuo, < t < T. (SI) 
Thus, the choice of the initial condition uq will be crucial in determining their impact. 

Now, it would be possible, in principle, to use the node of Juno to measure the 
gravito- magnetic effect. Indeed, the Lense-Thirring node precession is independent of 
i, while all the zonal precessions of Q vanish for i = 90 deg. It does not occur for the 
periJove and the mean anomaly, but they are not affected by the gravito-magnetic force 
for i = 90 deg. In reality, the situation will be likely different because of the unavoidable 
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orbit injection errors which will induce some departures Si of the Juno's orbital plane from 
the ideal polar configuration. Thus, unwanted, corrupting node zonal secular precessions 
will appear; their mis-modeling due to the uncertainties 6J£ may swamp the recovery of 
the Lense-Thirring effect if their determination by Juno will not be accurate enough. Note 
that there is no risk of some sort of a-priori 'imprint' effect of the Lense-Thirring effect 
itself on the values of the zonals retrieved from the Juno's periJove motion because the 
gravito-magnetic pericenter precession vanishes for polar orbits. 

By assuming the values quoted in Table [2] for the uncertainties 6Je, £ = 2,3, 4, 6, let us 
see what is the impact of a non-perfectly polar orbital geometry on the node Lense-Thirring 
precessions. The results are depicted in Figure [1] for each degree i separately; the initial 
condition uq = 90 deg has been used. It should be noted that, in view of the likely 
correlations among the determined zonals, a realistic upper bound of the total bias due to 
them can be computed by taking the linear sum of each mis-modeled terms. The major 



source of bias is the so far poorly known Jq, an improvement of four orders o 



magnitude. 



20041 ). would 



which sounds rather unlikely to be obtained even with Juno (lAnderson et al.l 
be required to push its aliasing effect at a percent level of the Lense-Thirring effect. The 
situation for the other zonals is more favorable; J4 should be known better than now by a 



Anderson et al 



fl2004h . Thus, 



factor 1000, which is, instead, a realistic goal according to 
we conclude that a nearly-polar orbit 1 deg off the ideal 90 deg case would likely prevent to 
obtain a measurement of the gravito-magnetic node precession at a decent level of accuracy. 

Thanks to the high eccentricity of the Juno's orbit, also the periJove and the mean 
anomaly are well defined, so that they can be used in a suitable way to remove the bias of 
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Node Q. only (co^= 90 deg) 



zi. 

to 



11000 
10500 
10000 
9500 
9000 
8500 
8000 
7500 
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6500 
6000 
5500 
5000 
4500 
4000 
3500 
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2500 
2000 
1500 
1000 
500 



- - J„ 



39 89.1 89.2 89.3 89.4 89.5 89.6 89.7 89.8 89.9 90 90.190.2 90.3 90.4 90.5 90.6 90.7 90.8 90.9 91 

/ (deg) 



Fig. 1. — Systematic percent bias on the Lense-Thirring node precession induced by the mis-modchng in 
the zonals J2, J3, J4, Jq according to Table [2]for 89 deg < i < 91 deg and ujq ~ 90 deg. 
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Jq and J2. Let us write down 

5tl = n,25J2 + (lJJ6 + finLT + ^n, (52) 



Su = U2SJ2 + UJ.6SJ6 + fJ'^LT + ^uj, (53) 



SM = M,2SJ2+ M,eSJe + Am- (54) 

Here denotes some sort of Observed- minus- Calculated (O — C) quantity for the 
rate of the Keplerian element \1/ which accounts for every unmo deled/mis-mo deled 
dynamical effects; it may be, for example, a correction to the modeled precessions to be 
phenomenologically estimated as a solve-for parameter of a global fit of the Juno's data. 



or it could be a computed time-series of 



"residuals" of by suitably overlapping orbital 



arcs. The gravito-magnetic force s 



lould be purposely not modeled in order to be fully 



present in 5\E'. The parameter fi 



ifll i 



in GTR and in Newtonian mechanics and accounts 
for the Lense-Thirring effect. The A terms include all the other systematic errors like the 
precessions induced by the mis-modeled parts of the second even zonal harmonic 6J4 and 
the first odd zonal harmonic 6Js, the mis-modeling due to the uncertainty in Jupiter's GM, 
etc. By solving for /i one obtains 

+ Ci5uj + C26M = tluT + Ci^LT + A, (55) 



^Since the Keplerian elements are not directly measurable quantities, we use here the 
term "residual" in an improper sense. 

^It is not one of the standard PPN parameters, but it can be expressed in terms of 7 as 
/x=(l + 7)/2. 
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with 



Cl 



M,6 - ^.6 M,2 
^.6 M.2 - M.6 ^.2 ' 



(56) 



C2 



(57) 



M.2 - M.Q UJ.2 

The combination of eq. (155|) . with eq. fl56|l and eq. f l57|) . is designed, by construction, to 
single out the combined Lense-Thirring precessions and to cancel out the combined seculac 
precessions due to J2 and Jq along with their mis-modeling. Instead, it is affected by A 
which acts as a systematic bias on the Lense-Thirring signal of interest. A globally includes 
the mis- modeled part of the combined precessions induced by J3 and J4; the sources of 
uncertainty reside in J3 and J4 themselves and in the Jupiter's GM through the mean 
motion n which enters the coefficients Clj,uj.i, Ai.i. 

In Figure [2] the impact of the mis-modeling in J3 and J4 for = 90 deg is depicted. 
In Figure [3] we use Uq = deg. In this case the situation is much more favorable because 
for a total departure of ±1 deg from i = 90 deg, an improvement of only two orders of 
magnitude in J3, which is, today, still compatible with zero, and J4 would be needed to 



reach the per cent level; 



the results by 



Jacobson 



et us recall that the expected improve ment in J4 with 



20031 ) is of three orders of magnitude (lAnderson et al.l 120041 ). Note 



resp ect to 



that a value of far from 90 deg is preferable to minimize the perturbations. 

Another potential source of systematic error is the Jupiter's GM whose uncertainty 
S{GM) indirectly affects eq. (155!) through the Keplerian mean motion n entering the 
uncanceled J3 and J4 combined precessions; 6n is also present via the mean anomaly itself. 
However, it turns out that it is of no concern because, according to the present-day level of 



^We include in them also the long-period harmonic terms for the reasons explained before. 
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Q. + c^co + c^M ((B^= 90 deg) 
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Fig. 2. — Systematic percent bias on the Lense-Thirring precessions, combined according to eq. ()55p . 
induced by the mis-modehng in the uncanceled zonals J3, J4 (Table [5]) for 89 deg < i < 91 deg and cjq = 90 
deg. 
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Fig. 3. — Systematic percent bias on the Lense-Thirring precessions, combined according to eq. ()55p . 
induced by the mis-modehng in the uncanceled zonals J3, J4 (Table [5]) for 89 deg < i < 91 deg and cjq = 
deg. 
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relative uncertainty f lJacobsoru |2003| ) 



its impact on the combined Lense-Thirring precessions is well below the percent level. 

2.2. A numerical approach 

Also in this case, we followed an alternative approach based on preliminary numerical 
simulations. We investigated the impact of the uncertainties in the first two jovian even 
zonals on a Juno's single six-hours pass by numerically simulating the probe's Doppler 
range-rate signals due to 5J2 and 5J4. By assuming for them values as large as 2 x 10~^° and 
3 X 10"^°, respectively, it turns out that the maximum Doppler shifts are roughly 1 — 1.5 
yUm s~^. Moreover, and more importantly, the time-dependent patterns of the even zonals' 
Doppler signals are quite different from the Lense-Thirring one removing the risk of an 
insidious mimicking bias. Another encouraging fact is that such simulations indicate that 
an inclination of even 91 deg would not compromise the recovery of the gravitomagnetic 
signal of interest. 



3. Discussion and conclusions 

The node VL of Juno, a recently approved spacecraft aimed to orbit Jupiter along a 
highly eccentric (rmin = 1.06-Rjup, rmax = 39-Rjup), polar {i = 90 deg) trajectory during one 
year to accurately map, among other things, its gravitational field, will be displaced by the 
general relativistic gravito-magnetic Lense-Thirring effect by about 572 m over the entire 
duration of the mission. 

We, first, explored the possibility of a high accuracy measurement of such an effect 
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by performing analytical calculations and interpreting them in a rather conservative 
fashion. A meter-level accuracy in determining the jovicentric orbit of Juno should not be 
an unrealistic goal to be reached. Equivalently, the gravito-magnetic node precession of 
Juno amounts to 68.5 mas yr~^, while the accuracy in measuring its node and peri Jove 
precessions should be of the order of 0.5 — 1 mas yr"'^, given the expected improvements in 
our knowledge of the departure of the jovian gravitational field from spherical symmetry. 
If the Juno's orbit was perfectly polar, the long-period node precessions induced by the 
zonal harmonics J^, £ = 2, 3, 4, 6, ... of the non-spherical jovian gravitational potential would 
vanish, thus removing a major source of systematic alias on the Lense-Thirring secular 
precession. In reality, unavoidable orbit injection errors will displace the orbital plane of 
Juno from the ideal polar geometry; as a consequence, the mis-modeled part of the node 
zonal precessions would overwhelm the relativistic signal for just 5i = ±1 deg, in spite of 
the expected improvements in = 2,4,6 of three orders of magnitude. A suitable linear 
combination of the node, the periJove u and the mean anomaly Ai will allow to cancel out 
the effect of J2 and Jq, the remaining uncanceled J3 and J4 will have an impact on the 
Lense-Thirring combined precessions which should be reduced down to the percent level or 
better by the improved low-degree zonals. 

Instead of looking at the cumulative, secular effects over the entire duration of the 
mission, we also followed an alternative approach by looking at single Doppler range-rate 
measurements over time spans six hours long centered on the the probe's closest approaches 
to Jupiter; it turned out that, in this way, the perspectives are even more favorable. We 
numerically simulated the characteristic Lense-Thirring pattern for a single science pass by 
finding a maximum value of the order of hundreds fim s~^, while the expected precision 
level in Juno's Doppler measurements is of the order of one /im s~^. Thus, by exploiting 
about 25 of the planned 33 total passes of the mission it would be possible to reach a 
measurement accuracy below the percent level. We repeated our numerical analysis also 
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for the Doppler range-rate signals of J2 and J4 by finding quite different patterns witli 
respect to the gravito-magnetic one; moreover, for a level of mismodeling of the order of 
2 — 3 X 10^^*^ in such zonals the maximum value of their biasing Doppler signals is about 
1 — 1.5 /im s"^. Our numerical analysis also shows that a departure from the nominal polar 
orbital geometry as large as 1 deg would not compromise the successful outcome of the 
measurement of interest, contrary to the conservative conclusions of our analytical analysis. 
Thus, this approach shows that there is not a high correlation between the Lense-Thirring 
parameter and the jovian gravity field parameters, although a covariance analysis would be 
needed to prove it. However, such a covariance analysis is outside the scope of the present 
paper. 

In conclusion, the potential error in the proposed Juno Lense-Thirring measurement is 
between 0.2 and 5 percent. Conversely, if one assumes the existence of gravito-magnetism 
as predicted by general relativity, the proposed measurement can also be considered as a 
direct, dynamical determination of the jovian proper angular momentum S by means of the 
Lense-Thirring effect at the percent level. 
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